Quantum vacuum properties of the intersubband cavity polariton field 
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We present a quantum description of a planar microcavity photon mode strongly coupled to a 
semiconductor intersubband transition in presence of a two-dimensional electron gas. We show that, 
in this kind of system, the vacuum Rabi frequency Qr can be a significant fraction of the inter- 
subband transition frequency W12. This regime of ultra-strong light-matter coupling is enhanced 
for long wavelength transitions, because for a given doping density, effective mass and number of 
quantum wells, the ratio Q.r/lu\2 increases as the square root of the intersubband emission wave- 
length. We characterize the quantum properties of the ground state (a two- mode squeezed vacuum) , 
which can be tuned in-situ by changing the value of Hr, e.g., through an electrostatic gate. We 
finally point out how the tunability of the polariton quantum vacuum can be exploited to generate 
correlated photon pairs out of the vacuum via quantum electrodynamics phenomena reminiscent of 
the dynamical Casimir effect. 



In the last decade, the study of intersubband elec- 
tronic transitions^ in semiconductor quantum wells has 
enjoyed a considerable success, leading to remarkable 
opto-electronic devices such as the quantum cascade 
lasersSiiM. In contrast to the more conventional inter- 
band transitions between conduction and valence bands, 
the frequency of intersubband transitions is not deter- 
mined by the energy gap of the semiconductor material 
system used, but rather can be chosen via the thickness of 
the quantum wells in the active region, providing tunable 
sources emitting in the mid and far infrared. 

One of the most fascinating aspects of light-matter 
interaction is the so-called strong light-matter coupling 
regime, which is achieved when a cavity mode is reso- 
nant with an electronic transition of frequency u>\2, and 
the so-called vacuum Rabi frequency exceeds the cav- 
ity mode and electronic transition linewidths. The strong 
coupling regime has been first observed in the late '80s 
using atoms in metallic cavities^, and a few years later 
in solid-state systems using excitonic transitions in quan- 
tum wells embedded in semiconductor microcavitiesi In 
this regime, the normal modes of the system consist of 
linear superpositions of electronic and photonic excita- 
tions, which, in the case of semiconductor materials, are 
the so-called polaritons. In both these systems, the vac- 
uum Rabi frequency flu does not exceed a very small 
fraction of the transition frequency u>i2- 

Recently, Dini et alM- have reported the first demon- 
stration of strong coupling regime between a cavity pho- 
ton mode and a mid-infrared intersubband transition, in 
agreement with earlier semiclassical theoretical predic- 
tions by LiuS. The dielectric Fabry-Perot structure real- 
ized by Dini et al& consists of a modulation doped mul- 
tiple quantum well structure embedded in a microcavity, 
whose mirrors work thanks to the principle of total in- 
ternal reflection. The strong coupling regime has been 
also observed in quantum well infra-red detectors^. As 
we will show in detail, an important advantage of using 
intersubband transitions is the possibility of exploring 
a regime where the normal-mode polariton splitting is 



a significant fraction of the intersubband transition (in 
the pioneering experiments by Dini et al&, 2hfl]i = 14 
meV compared to huji2 — 140 meV). Furthermore, re- 
cent experiments have also demonstrated the possibility 
of a dramatic tuning of the strong light-matter coupling 
through application of a gate voltag ell which is able to 
deplete the density of the two-dimensional electron gas. 

Although the quest for quantum optical squeezing ef- 
fects in the emission from atoms strongly coupled to a 
cavity mode has been an active field of research^, all sys- 
tems realized up to now show a vacuum Rabi frequency 
flu much smaller than the frequency of the optical tran- 
sition. In this parameter regime, the relative importance 
of the anti-resonant terms in the light-matter coupling 
is small and, as far as no strong driving field is present, 
they can be safely neglected under the so-called rotating- 
wave approximation. In the presence of a strong driving 
field, however, anti-resonant terms are known to play a 
significant role, giving, e.g., the so-called Bloch-Siegert 
shift in magnetic resonance experiments^, or determin- 
ing the quantum statistical properties of the emission 
from dressed-state lasers^. 

A few theoretical studies have pointed out the intrin- 
sic non-classical properties of exciton-polaritons in solid- 
state systemsifiiiSiii, but the small value of the ratio 
Q,n/ui exci typically less than 0.01, has so far prevented the 
observation of quantum effects due to the anti-resonant 
terms of the light-matter coupling. All the squeezing 
experiments that have been performed so far in fact re- 
quired the presence of a strong coherent optical pump 
beam in order to inject polaritons and take advantage of 
nonlinear polariton parametric processesi2ii2i22i2ii22. 

In this paper, we show that in the case of intersub- 
band cavity polaritons, it is instead possible to achieve 
an unprecedented ultra-strong coupling regime, in which 
the vacuum Rabi frequency f2fl is a large fraction of the 
intersubband transition frequency uj\2- To this purpose, 
transitions in the far infrared are most favorable, because 
the ratio SIr/w^ scales as the square root of the inter- 
subband emission wavelength. Within a second quanti- 
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zation formalism, we characterize the polaritonic normal 
modes of the system in the weak excitation limit, in which 
the density of intersubband excitations is much smaller 
than the density of the two-dimensional electron gas in 
each quantum well (in this very dilute limit, the inter- 
subband excitations behave as bosons). We point out 
the non-classical properties of the ground state, which 
consists of a two-mode squeezed vacuum. As its proper- 
ties can be modulated by applying an external electro- 
static bias, we suggest the possibility of observing quan- 
tum electrodynamics effects, such as the generation of 
correlated photon pairs from the initial vacuum state. 
Such an effect closely reminds the so-called dynamical 
Casimir effectSiSiiS^, whose observation is still an open 
challenge and is actually the subject of intense effort. 
Many theoretical works have in fact predicted the gener- 
ation of photons in an optical cavity when its properties, 
e.g. the length or the dielectric permittivity of the cavity 
spacer material, are modulated in a rapid, non-adiabatic 
way^LSS. 

The present paper is organized as follows. In Sec. Hwc 
describe the system under examination and in Sec. ITllwe 
introduce its Hamiltonian. The scaling of the coupling in- 
tensity with the material parameters is discussed in Sec. 
IIIII while Sec. II VI is devoted to the diagonalization of the 
Hamiltonian and the discussion of the polaritonic nor- 
mal modes of the system in the different regimes. The 
quantum ground state is characterized in Sec. [3 and 
its quantum properties are pointed out. Two possible 
schemes for the generation of photon pairs from the ini- 
tial vacuum by modulating the properties of the ground 
state are suggested in Sec. IVII Conclusions are finally 
drawn in Sec. IVIII 




FIG. 1: (a) Sketch of the considered planar cavity geom- 
etry, whose growth direction is called z. The cavity spacer 
of thickness L ca , v embeds a sequence of uqw identical quan- 
tum wells. The energy of the cavity mode depends on the 
cavity photon propagation angle 0. (b) Each quantum well 
contains a two-dimensional electron gas in the lowest subband 
(obtained through doping or electrical injection). The transi- 
tion energy between the first two subbands is huii2. Only the 
TM-polarized photon mode is coupled to the intersubband 
transition and a finite angle 6 is mandatory to have a finite 
dipole coupling, (c) Sketch of the energy dispersion Ei(q) and 
Evil) — Ei(q) + hu>i2 of the first two subbands as a function 
of the in-plane wavevector q. The dispersion of the inter- 
subband transition is negligible as compared to the one of the 
cavity mode. For a typical value of the cavity photon in-plane 
wavevector k, one has in fact E2(\k + q\) — Ei(q) ~ hui2- 



I. DESCRIPTION OF THE SYSTEM 

In the following, we will consider a planar Fabry-Perot 
resonator embedding a sequence of uqw identical quan- 
tum wells (see the sketch in Fig. QJ. Each quantum well 
is assumed to be doped with a two-dimensional density 
N2DEG of electrons, which, at low temperatures, popu- 
late the first quantum well subband. Due to the presence 
of the two-dimensional electron gas, it is possible to have 
transitions from the first to the second subband of the 
quantum well. We will call 7kJi2 the considered intersub- 
band transition energy. If we denote with z the growth 
direction of the multiple quantum well structure, then 
the dipole moment of the transition is aligned along z, 
i.e., di2 = di2Z. This property imposes the well known 
polarization selection rule for intersubband transitions in 
quantum wells, i.e., the electric field must have a com- 
ponent along the growth direction. We point out that in 
the case of a perfect planar structure, the in-plane wave- 
vector is a conserved quantity, unlike the wave-vector 
component along the z direction. Therefore, all wave- 
vectors k will be meant as in-plane wave-vectors, unless 
differently stated. 



In the following, we will consider the fundamental cav- 
ity photon mode, whose frequency dispersion is given by 
<^cav,fc = ~^=y/kz + where is the dielectric con- 
stant of the cavity spacer and k z is the quantized photon 
wavevector along the growth direction, which depends 
on the boundary conditions imposed by the specific mir- 
ror structures. In the simplest case of metallic mirrors, 
k z — j^—, with i cav the cavity thickness. 



II. SECOND QUANTIZATION HAMILTONIAN 

In this Section, we introduce the system Hamiltonian 
in a second quantization formalism. In the following, 
we will call a k the creation operator of the fundamental 
cavity photon mode with in-plane wave- vector k. Note 
that, in order to simplify the notation, we will omit the 
polarization index of the photon mode, which is meant 
to be Transverse Magnetic (TM)-polarized (also known 
as p-polarization). This photon polarization is neces- 
sary to have a finite value of the electric field compo- 
nent along the growth direction z of the multiple quan- 
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turn well structure, direction along which the transition 
dipole of the intersubband transition is aligned, b k will 
be instead the creation operator of the bright intersub- 
band excitation mode of the doped multiple quantum 
well structure. In the simplified case of uqw identical 
quantum wells that are identically coupled to the cavity 
photon mode, the only bright intersubband excitation is 
the totally symmetric one, with an oscillator strength 
tiqw times larger than the one of a single quantum well. 
The uqw — 1 orthogonal excitations are instead dark and 
will be neglected in the following. The creation operator 
corresponding to the bright intersubband transition can 
be written as 
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where N%d eg is the density of the two-dimensional elec- 
tron gas in each quantum well and S is the sample area. 
The fermionic operator cf q annihilates an electron be- 
longing to the first subband and j-th quantum well, 
while C2 g +k creates an electron in the second subband 
of the same well, fc^ is the Fermi wavevector of the two- 
dimensional electron gas, whose electronic ground state 
at low temperature is 



w=n n 
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3=1 \q\<k F 



where \0} CO nd is the empty conduction band state. 

In the following, we will consider the situation of a 
weakly excited intersubband transition, i.e., 



i E< 6 ^> « N ' 
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2DEG 



(3) 



In this dilute limit, the intersubband excitation field is 
approximately bosonic, namely 



[bk,b 



8k,k' 



(4) 



Starting from the coupled light-matter Hamiltonian 
of the semiconductor and retaining only the consid- 
ered cavity photon mode for the electromagnetic field 
and the considered intersubband transition for the elec- 
tronic polarization field, one finds a standard Hopficld- 
likc HamiltonianSS 



H — Hq + H r 



Ha 



(5) 



which consists of three qualitatively different contribu- 
tions, namely 
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H ° = E ^ w cav,fc f a> k a k + ^ j + E huJl2 b k bk ' ( 6 ) 

k ^ ' k 



Hres = ^E ( a fc^ _ a * b l) 

k 

+D k (a k a k + afeajfe) | , (7) 



HanU = ^E {^R'k _ a k b -k) 

k 

+D k (a k a_ k + a^,al fc ^| . (8) 

Hq in Eq. |JBJ describes the energy of the bare cavity 
photon and intersubband polarization fields, which de- 
pend on the numbers a k a k , b^ k b k of cavity photons and 
intersubband excitations, respectively. 

H res in Eq. JJJ) is the resonant part of the light-matter 
interaction, depending on the vacuum Rabi energy hQRk 
and on the related coupling constant Dk ■ The terms pro- 
portional to describe the creation (annihilation) of 
one photon and the annihilation (creation) of an inter- 
subband excitation with the same in-plane wavevector. 
In contrast, the term proportional to Dk contains only 
photon operators, because it originates from the squared 
electromagnetic vector potential part of the light- matter 
interaction. Note that this term in H res depends on the 
photon number operator a k a k as the bare cavity pho- 
ton term in Eq. 10. Hence, it gives a mere blueshift 
(Dfc > 0) of the bare cavity photon energy 7La> C av,fc- 

Finally, H ant i in Eq. JSJ contains the usually neglected 
anti-resonant terms, which correspond to the simultane- 
ous destruction or creation of two excitations with op- 
posite in-plane wavevectors. The terms proportional to 
f2.H,fc describe the creation (or destruction) of a cavity 
photon and an intersubband excitation, while the terms 
proportional to D k describe the corresponding process 
involving a pair of cavity photons. 

Before continuing our treatment, we wish to point out 
that the considered Hamiltonian in Eq.® contains only 
the energy associated to the fundamental cavity mode 
(including the zero-point energy ^ fc | fi.w C av,fc), the en- 
ergy associated to the creation of intersubband excita- 
tions and the full light-matter interaction between the 
considered modes. The energy terms associated to the 
other photon modes, the electronic energy of the filled 
electronic bands as well as the electrostatic energy as- 
sociated to an applied bias have been here omitted for 
simplicity, as they do not take part in the dynamics dis- 
cussed in the following of the paper. 

III. SCALING OF THE INTERACTION 

The specific values of the coupling constants Qr^ and 
Dk depend on the microscopic parameters of the inter- 
subband microcavity system. 

The so-called vacuum Rabi energy UQr^ is the Rabi 
energy obtained with the electric field corresponding to 
one photon^S. For the system under consideration^, 
the polariton coupling frequency for the TM-polarized 
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where is the dielectric constant of the cavity, L°f v the 
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FIG. 2: Coupling ratio as a function of the intersub- 

band emission wavelength A12 (/im). Parameters: /12 = 12.9 
(GaAs quantum well), cavity spacer refraction index yfe^, — 
3, n e Q W = 50, N 2 deg = 5 x 10 11 cm" 2 and = 60°. The 
Fabry-Perot resonator is a A/2-microcavity. Results obtained 
from the analytical expressions in Eqs. 1141 and 1151 . 



effective thickness of the cavity photon mode (which de- 
pends non-trivially on the boundary conditions imposed 
by the specific mirror structures), and n C Q W the effective 



number of quantum wells (^? ff 



in the case of 



"QW - n QW 

quantum wells which are identically coupled to the cav- 
ity photon field and which are located at the antinodes 
of the cavity mode electric field) . The oscillator strength 
of the considered intersubband transition reads 



/12 = 2m Q u 12 d\ 2 /'h 



(10) 



where too is the free electron mass and d\i is the electric 
dipole moment of the transition. Under the approxima- 
tion of a parabolic energy dispersion of the quantum well 
subbands, the oscillator strengths of the different inter- 
subband transitions satisfy the /-sum rula^ 
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flj = m /m* 



(11) 



where to* is the effective electron mass of the conduction 
band. In particular, for our case of a deep rectangular 
well, the sum rule is almost saturated by the first inter- 
subband transition /12 ~ mo/m*. Finally, 9 is the prop- 
agation angle inside the cavity (which is different from 
the propagation angle in the substrate) , and is related to 
the in-plane wavevector k by k/k z = sinf?/ cos 9. 

As we will see in the next section, the relevant pa- 
rameter quantifying the importance of the quantum ef- 
fects considered in this paper is the dimensionless ratio 
^fl,fc I . e s/ a; i2j where k res is the resonance in-plane wavevec- 
tor such as fi.Wcav,fc ros = naji2- In the system studied 
by Dini et al&, this ratio is already significant, namely 
&R,k/ui2 — 0.05. Here, we show that the ratio FIr^/wiv 
can be largely increased designing structures in the far 



infra-red, by increasing the number of quantum wells and 
by choosing semiconductors with smaller effective mass. 

Let be 8 ICS the cavity propagation angle corresponding 
to fc res - From the relation 



, sin 9 r . 



we get that for metallic mirrors 

A12 



L, 



,cos 9 r 



(12) 



(13) 



where 2n/Xi 2 = w 12 /c. Under these conditions, the light 
matter coupling ratio at the resonance angle is 

&R,k 



W12 



with 



'/ : 



/12 sin 2 9 ICS cos fl rcs N 2DEG n QW 



(14) 



(15) 



Note that the prefactor given in Eq. l|15|) has a weak de- 
pendence on A12. In fact, in the limit case of a rectan- 
gular quantum well with high potential barriers, fi 2 = 
0.96 mo/m* and does not depend at all on A12. More 
refined calculations^ including the non-parabolicity of 
the semiconductor band and the finite depth of the po- 
tential well show that /12 has a moderate dependence on 
the emission wavelength A12 (it actually increases with 
A12). Hence, the normalized vacuum Rabi frequency 
f2fl,fc re s/ a; i2 increases at least as \/Ai2- The predictions 
of Eqs. and (|T5|) are reported in Fig. [21 for a sys- 

tem of 50 GaAs quantum wells and a doping density 
N 2 deg — 5 x 10 11 cm" 2 . For an intersubband emission 
wavelength of lOO^rn, the ratio f2#/a->i2 can be as high as 
0.2. The values in Fig. [21 can be significantly increased 
using semiconductors with smaller effective mass, such as 
InGaAs/ AlInAs-on-InB2i. 

To complete our description, we need to provide the ex- 
plicit expression for the coefficient D^, which quantifies 
the effect of the squared electromagnetic vector potential 
in the light-matter interaction. Generalizing Hopfield's 
procedure 2 ^ to the case of intersubband transitions, we 
find that all the intersubband transitions give a contri- 
bution to Dk, namely 



D k = 



T,j hi ^R,k 



(16) 



/12 ^12 

However, as the oscillator strength of a deep rectangu- 
lar well is concentrated in the lowest transition at W12, 
the effect of the higher transitions is a minor correction, 
namely 



D u 



1.04 



O 2 



o 2 



(17) 



UJi 2 UJi 2 

Note that for a quantum well with a parabolic confine- 
ment potential V(z) — (l/2)m*uj 2 2 z 2 , the expression 
Dk = fl 2 R k /oJi2 would be exact, since in this case all the 
intersubband oscillator strength is exactly concentrated 
in the lowest transition ui 2 . 



5 



2.5 




FIG. 3: Normalized polariton frequencies uoLP,k/ 'wi2 and 
Luup,k/ijJi2 as a function of 0,R,k/ui2 for Dk = fl% >k /iJi2- The 
calculation has been performed with tj cav ,fc = u>i2- Note that 
for a given microcavity system, Q.R,k/u\2 can be tuned in-situ 
by an electrostatic bias, which is able to change the density 
of the two-dimensional electron gas. 



IV. INTERSUBBAND POLARITONS 

As all the terms in the Hamiltonian H — Hq + H res + 
Hanti are bilinear in the field operators, H can be exactly 
diagonalized through a Bogoliubov transformation. Fol- 
lowing the pioneering work by Hopfield^, we introduce 
the Lower Polariton (LP) and Upper Polariton (UP) an- 
nihilation operators 

PjM = Wj,k afc + %j,k b k + Vj,k a)_ k + z j>k fri fc , (18) 

where j E {LP, UP}. The Hamiltonian of the system 
can be cast in the diagonal form 

H = E G + ]T X l> : ,Kl'> « ■ ( 19 ) 

j£{LP,UP} k 

where the constant term Eg will be given explicitly later. 
The Hamiltonian form in Eq. (|19fl is obtained, provided 
that the vectors 

v hk = (w hk ,x-j.k,yj.k,Zj^) T (20) 

satisfy the eigenvalues equation 

M k Vj, k = uJj,kVj,k (21) 

with ujj^ > 0. The Bose commutation rule 

[Pi,*»Pi',fe'] = § jj' 6 k,k' (22) 
imposes the normalization condition 

w lk w j',k + x*j,k x j',k ~ vlkVrM - Zj,k z j',k = Sj,j' ■ (23) 




FIG. 4: Mixing fractions for the Lower Polariton (LP) 
mode as a function of Q.R,k/u\2 (see Eq. HSU in the text). 
The calculation has been performed for the resonant case 
Wcsv.t = 1^12 , as in the previous figure. Panel (a): \wLP,k\ 2 
(thin solid line), |a^z,p,fe| 2 (thick solid line). Note that for 
O.R,k/Lo 12 < 1, \w LP m\ 2 ~ \x LP ,k\ 2 ~ 1/2. Panel (b): 
\yLP,k\ 2 (thin dashed line), \zLP,k\ 2 (thick dashed line). For 
< 1, \VLP,k\ 2 ^ \z L p,k\ 2 0. The Upper Polariton 
(UP) fractions (not shown) are simply \wup.k\ 2 = \xLP.k\ 2 , 

1 2 ! 1 2 I 1 2 i |2" i |2 1 1 2 

\XUP,k\ = \W L P,k\ , \VUP,k\ = \ZLP,k\ , \ZuP,k\ = \VLP,k\ ■ 



The Hopfield-like matrix for our system reads 



M k = 



1 






-2D k 


-i^R,k 




i^R,k 


W12 


-i^R,k 







2D k 


-i^R,k 


—^cav.k — 2-Dfe 


— i^R,k 


\ 







i^R,k 


-Wl2 



The four eigenvalues of M k arc {±07LP,fe, ±<^t7P,fc}- Under 
the approximation D k = Vl 2 R k /uj\2 (i.e., all the oscillator 
strength concentrated on the U)t2 transition), detMfc = 
(^cav.fe W12) 2 , giving the simple relation 

L^LP,k LdUP,k — ^12 ^cav,/c 5 

(25) 

i.e., the geometric mean of the energies of the two po- 
lariton branches is equal to the geometric mean of the 
bare intersubband and cavity mode energies. The de- 
pendence of the exact polariton eigenvalues as a function 
of f2ii,fc/wi2 is reported in Fig. [31 for the resonant case 



A. Ordinary properties in the limit Q,R lk / '^12 <C 1 

In the standard case ^lp n k/ ^12 <C 1, the polariton op- 
erator can be approximated as 

pj.k ~ uij t k a k + Xj >k b k , (26) 
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with |u>j,fc| 2 + |a;j,fc| 2 — 1- This means that the annihi- 
lation operator for a polariton mode with in-plane wave- 
vector k is given by a linear superposition of the photon 
and intersubband excitation annihilation operators with 
the same in-plane wavevector, while mixing with the cre- 
ation operators (represented by the coefficients and 
Zj_k) is instead negligible [see Fig. |3J. In this limit, the 
geometric mean can be approximated by the arithmetic 
mean and Eq. (|25|l can be written in the more usual 
form: 

^LP,k + &UP,k — ^cav.fc + ^>12 ■ (27) 

For the specific resonant wavevector fc rcs such that 
^cav.fc rca = W12, the polariton eigenvalues are 

w LP(!7P),fc rra — ^12 T ^fi.fcrcs ' (28) 

and the mixing fractions are \wLPk , | 2 — \%LPk | — 
1/2- 

B. Ultra-strong coupling regime 

When the ratio ^Ir^/oju is not negligible compared to 
1, then the anomalous features due to the anti-resonant 
terms of the light-matter coupling becomes truly rele- 
vant. 

In the resonant w cav ,j ![os = UJ12 case and under the ap- 
proximation Dk — fl 2 Rk /uJi2, the polariton frequencies 
are given by 



ULP(UP),k IOB = y W12 + {^R,k Tea ) 2 T ^R,k tea , (29) 

which, as it is apparent in Fig. [3J corresponds to 
a strongly asymmetric anti-crossing as a function of 
^i?,fc rca / w i2- This is due to the combined effect of the 
blue-shift of the cavity mode frequency due to the terms 
proportional to D k in Eq. (jTJ, and of the anomalous 
coupling terms in Eq.©. 

These same effects contribute to the non-trivial evo- 
lution of the Hopfield coefficients shown in Fig. 0] The 
anomalous Hopfield fractions \yLP,k\ 2 and \zLP,k\ 2 signif- 
icantly increase because of the anomalous coupling, and 
eventually become of the same order as the normal ones 
|^LP.fc| 2 and \wLP,k\ 2 - Due to the normalization condi- 
tion 

K fe | 2 + |^ fe | 2 -|%,fcl 2 -M 2 = i , (30) 



this affects the ordinary fractions \wLP.k\ 2 , \xLP,k\ 2 as 
well. Owing to the blue-shift of the cavity photon fre- 
quency induced by the light-matter coupling, at the 
resonance wavevector k — k rcs the lower polariton be- 
comes more matter-like (i.e., \xLP,k ros \ 2 > \wLP,k,c S \ 2 and 
\zi,p,k IOS \ 2 > \yLP,k IOS | 2 ), while the upper polariton more 
photon-like. In this resonant case, the UP Hopfield coef- 
ficients (not shown) are simply related to the LP ones 
by: \w UP m icb \ 2 = \x L p,k rcB \ 2 , \xup,k IO J 2 = |wLRfc ro J 2 , 
\yup,k r J 2 = \zLP,k,J 2 , \zup,k r J 2 = \VLP,k,J 2 - 

V. THE QUANTUM GROUND STATE 

A. The normal vacuum state |0) for = 

In the case fi^ = (negligible light-matter interac- 
tion), the quantum ground state \G) of the considered 
system is the ordinary vacuum |0) for the cavity photon 
and intersubband excitation fields. Such ordinary vac- 
uum satisfies the relation 

o fc |0) = 6 fc |0) -0 , (31) 

which means a vanishing number of photons and inter- 
subband excitations: 

(0|4a fe |0) - (0|&t& fc |0) = <0|4& fe |0) = (32) 

and no anomalous correlations, i.e., 

(0\a k a k , |0) = (0\b k b k , |0) = (0|a fc fo fc ,|0> = . (33) 

B. The squeezed vacuum state 

With a finite tl^k, the ground state of the system \G) 
is no longer the ordinary vacuum |0) such that: 

o fc |0) - & fc |0) = , (34) 

but rather the vacuum of polariton excitations: 

Pj , k \G) = 0. (35) 

As the polariton annihilation operators are linear su- 
perpositions of annihilation and creation operators for 
the photon and the intersubband excitation modes, the 
ground state \G) is, in quantum optical terms, a squeezed 
stat o 34 ' 35 . By inverting Eq. (JTSJl , one gets 



/ a k \ 
b k 

V bl k J 



V- 



"LP,k 
c LP,k 
~VLP,k 
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W 
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XLP,k %UP,k ) 
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VhP-k 
V PlJP.-k I 
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a time-dependent external electrostatic bias. In partic- 
ular, we shall discuss how this remarkable tunability of 
the system can be used to "unbind" the virtual photons 
by modulating the parameters of the system in a time- 
dependent way, and generate some radiation which can 
be actually detected outside the cavity. These issues will 
be the subject of sec. I VII 



C. The ground state energy 

Also the energy Eg of the quantum ground state has a 
significant dependence on the coupling flR,k- Defining Eq 
as the ground state energy of the uncoupled (£lR, k = 0) 
system, we have that: 



FIG. 5: Solid line: normalized differential zero-point energy 
(per mode) Au>zp(k)/uJi2 as a function of Dashed 
line: Dk/uJ\2- Dotted line: normalized correlation contribu- 
tion Ao;|°p r (fc) /wi2. The calculation has been performed with 

I^cav.fc = Wl2- 



from which, using Eq. (|35() and the boson commutation 
rules, we obtain that the ground state contains a finite 
number (per mode) of cavity photons and intersubband 
excitations: 

(G\a{a k \G) = \y LP , k \ 2 + \y UP ^\ 2 (37) 
(G\b{b k \G) = \z LPtk \ 2 + \z UP , k \ 2 , (38) 

as well as some correlation between the photon and in- 
tersubband fields: 

(G\a< k b k \G) = ytp, k z L p,k + Vup,k z up,k ■ (39) 

Moreover, significant anomalous correlation exist be- 
tween opposite momentum components of the fields: 

(G\a k a- k \G) = -w* LP>k y L p,k - Wup,kVup,k (40) 
(G\b k b- k \G) = -x* LPtk z L p,k ~ x*up,k z up,k (41) 
(G\b k a- k \G) = -x* LPtk yLP,k ~ xhp,kVUP,k ■ (42) 

Note that the finite photonic population which is present 
in the ground state \G) of our system is composed of 
"virtual" photons. In the absence of any perturbation or 
modulation of the parameters of the quantum Hamilto- 
nian, these virtual photons can not escape from the cavity 
and therefore do not result in any observable emitted ra- 
diation (indeed, energy would not be conserved in such 
a process). 

As it has been shown in Fig. 0] the Hopfield coeffi- 
cients Xj^k, Vj,ki w j,k, Zj,k as well as the ground state \G) 
of the system strongly depend on the vacuum Rabi energy 
f2ii,fc, which in our case can be dramatically modulated 
in situ, e.g. by changing the electron density N2DEG via 



E G - E = [KD k - J2 h ^A\v 3 M 2 + K*| 2 )] ■ 

k je{LP,UP} 

(43) 

Note that this energy difference includes only the contri- 
bution of the zero-point fluctuations of the intersubband 
polariton field and does not take into account the other 
contributions coming, e.g. from the change of the elec- 
trostatic energy of the system (which is imposed by an 
applied bias), as already discussed at the end of Sec. ITTI 
The (always positive) term D k in Ea. (|43|l is the zero- 
point energy change due to the mere blue-shift of the 
bare cavity mode frequency and does not correspond to 
any squeezing effect. The second term is instead due to 
the mixing of creation and annihilation operators into 
the polaritonic operators as described in Ea. (|18f) and is 
proportional to the number of virtual photons and inter- 
subband excitations present in the ground state \G) of the 
system according to Eas. (|37l38ll . As it is usual for a cor- 
relation contribution, it tends to lower the ground state 
energy. It is interesting to study the differential "zero- 
point" energy per mode HAcuzp(k), whose sum over all 
the fc-modes gives the quantum ground state energy dif- 
ference Eg — Eq . The differential " zero-point" frequency 
reads 

Auj ZP (k)=D k + AuTf (k) , (44) 
where the (negative) correlation contribution reads 

Au J %T(k) = - £ Wjlfc (\ yj , k \ 2 + \z 3 M 2 ) ■ (45) 
je{LP,UP} 

These quantities (normalized to 0J12) are plotted in Fig. 
[3]as a function of ^l Ptk /uJi2 for the resonant case u> C av,fc = 
oj\2 ■ Although it is the diagonal blueshift which gives the 
dominant contribution to the ground state energy shift, 
the negative contribution due to the correlation effects 
is important, being as large as — 0.13?kJi2 already for 

MR,k/ui2 = 0.5. 
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X 

FIG. 6: Sketch of a possible set-up for the generation of cor- 
related photon pairs in the intersubband cavity system. The 
vacuum Rabi frequency of the intersubband cavity system can 
be modulated through an electric gate, which changes the den- 
sity of the two-dimensional electron gas or alternatively the 
dipole moment of the intersubband transition. A modulation 
of the bias is expected to induce the emission of correlated 
photon pairs with opposite in-plane wavevectors. This kind 
of radiation can be optimally guided out of the cavity through 
wedged lateral facets, with inclination equal to the resonance 
angle res . 



VI. TUNING THE QUANTUM VACUUM: 
QUANTUM RADIATION EFFECTS 

The possibility of tuning in a dramatic way the prop- 
erties (energy and squeezing) of the ground state of the 
system, as well as the significant number of (virtual) ex- 
citations already present in the ground state suggests 
that the present system could be a potential labora- 
tory to study Quantum Electro-Dynamics (QED) phe- 
nomena, which are reminiscent of the dynamic Casimir 
effect a 23 i 24 i 25 . In particular, we shall discuss how a time- 
modulation of the ground-state properties of the sys- 
tem can parametrically produce real excitations above 
the ground state of our cavity, which then escape from 
the cavity as photons and propagate in the external free- 
space. 

In the typical arrangement for the observation of the 
dynamical Casimir effect, one has to modulate in time 
the properties of an optical cavity and, in particular, its 
resonance frequencies. Several proposals have appeared 
in order to do this: in the simplest ones, one has to peri- 
odically move the mirrors so as to modify the boundary 
conditions of the field— i— Other proposals^ deal with 
a time-dependance of the refractive index of a dielectric 
medium placed inside the cavity. A recent work proposes 
to vary the effective length of the cavity by changing the 
reflectivity of a composite mirror—. 

The main peculiarity of our system as compared to 
previous proposals is due to the possibility of modulat- 
ing the properties of the ground state in a much stronger 
way, due to the ultra-strong and tunable light-matter 
coupling. 

In the next subsection, we shall give a detailed analy- 
sis of a simple gedanken experiment, where the vacuum 
Rabi frequency is assumed to be switched off in an in- 



stantaneous way. This scheme has the merit of allowing 
to grasp the essential physics of the problem, providing 
quantitative estimates without the need of embarking in 
complicate calculations. 

A complete and quantitative calculation of the spec- 
tral shape and intensity of the emitted radiation for the 
most relevant case of a periodic modulation of Qr^u is be- 
yond the scope of the present paper, as it would require 
a careful analysis of the coupling of the cavity system to 
the extra-cavity field as well as of the other non-radiative 
loss mechanisms of the electronic system— This is work 
actually in progress and here we shall restrict ourselves 
to a very qualitative discussion of its main features. 

A. Abrupt switch off of the vacuum Rabi energy 

Let us suppose that the considered intersubband cavity 
system is in the ground state |G). As we have already dis- 
cussed, the squeezed vacuum |G) contains a finite number 
of cavity photons and intersubband excitations because 
of the correlations due to the anomalous coupling terms 
in Eq.©. 

If one switches off the vacuum Rabi frequency fln.k 
of the system in an abrupt, non-adiabatic way by sud- 
denly depleting the electron gas, the photon mode does 
not have the time to respond to the perturbation and 
will remain in the same squeezed vacuum state as before. 
As this state is now an excited state of the Hamiltonian 
for fl^k = 0, the system will relax towards its ground 
state, which now corresponds to the standard vacuum, 
by emitting the extra photons as propagating radiation. 

One possible way to collect this quantum vacuum ra- 
diation is through the set-up sketched in Fig. which 
allows one to collect the photons which are emitted with 
internal propagation angle 9 around the resonance value 
Ores- If one neglects the losses due to the background 
absorption by the dielectric material forming the micro- 
cavity, an estimate of the number of emitted photons 
can be obtained as follows. The number of photon states 
(per unit area) in the 2D momentum volume d 2 k is sim- 
ply d 2 k/(2ir) 2 . Hence, the differential density of photons 
(per unit area) in the 2D momentum volume d 2 k is 

d 2 k 

dpphot = (G\a{a k \G) , (46) 

where the photon number (G|a^,afc|G) in the quantum 
ground state is given by Eq. (|37|l . Now, all the expecta- 
tion values depend only on |fc| and hence we can rewrite 
the momentum volume as d 2 k = 2irkdk. Knowing that 
the in-plane wavevector k is given by the relationship 
k = k z tan(6>) and using Eq. (|12|l . we find the final result 

^r(0 res ) = i^e 00 tan(e res )(G|aia fe |G) . (47) 
do 2n c l 

To give a numerical application of Eq. (|4TJl . let us 
consider an intersubband cavity system with hu>\2 = 140 
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meV, resonance angle 9 1CS = 65° and hflR t k lea = 7 meV 
(these are approximately the values in the sample mea- 
sured by Dini et aft). For these parameters, Eq. 147(1 
gives the differential photon density dp p hot/dO ~ 1 x 10 5 
cm~ 2 rad -1 . 

Note that the emission corresponding to the fc-modc 
is correlated to the emission corresponding to the mode 
with opposite in-plane wavevector, as shown in Eq. (|40|l . 
Indeed, the "quantum vacuum radiation" here described 
consists in the emission of correlated photon pairs2£. 

B. Periodic modulation of flR.k 

The requirement of a abrupt, non-adiabatic, switch-off 
of the Rabi coupling QR.k imposes very stringent limits 
on the time-scale t sw over which the electrostatic bias has 
to be applied. In particular, we expect that in order to 
maximize the quantum vacuum radiation generation, t sw 
can not be too much longer that the oscillation period of 
the lower polaritonic mode. 

It is then perhaps more accessible from an experimen- 
tal point of view to try to detect the vacuum radiation 
by periodically modulating the vacuum Rabi frequency 
at an angular frequency u> mo d 

Q,R,k{t) = &R,k + AQR,k sin(uj mod t). (48) 

Note that in principle this kind of modulation can be ob- 
tained not only through a gate-induced depletion of the 
two-dimensional electron gasii, but also by modulating 
the dipole moment of the intersubband transition or al- 
ternatively the reflectivity of the mirrors. As all the rele- 
vant physical quantities in the present problem (polariton 
energies, Hopfield coefficients, ground state energy) de- 
pend in a nonlinear way on the vacuum Rabi frequency 
^^R,k, we expect that for large modulation amplitudes 
high order harmonics of the fundamental modulation fre- 
quency LUmod will play a significant role in the paramet- 
ric process which is responsible for the vacuum radiation 
generation^. In particular, emission will be enhanced if 

Uj >h + Uj',-h = r uJ m od , (49) 

with r being a generic positive integer number, and 
j, f = {LP,UP}. This is the phase-matching condition 
for the parametric generation of two polaritons with op- 
posite momentum. As usual, the narrower the polaritonic 
resonance, the stronger the resonant enhancement. 

As it is generally the case for parametric processes in a 
cavity, the number of photons which are generated in the 
cavity and then emitted as radiation is determined by a 
dynamical equilibrium between the parametric processes 
generating them and the losses, the radiative as well as 



the non-radiative ones^ 6 .. For a complete and quantita- 
tive treatment of these issues, further investigations are 
in progress. 

VII. CONCLUSIONS 

In conclusion, we have shown that in the intersubband 
cavity polariton system, a new regime of ultra-strong cou- 
pling can be achieved, where the vacuum Rabi frequency 
Qr, is a large fraction of the intersubband transition fre- 
quency uji2 . This scenario appears to be easier to achieve 
in the far infrared, since the ratio Ur/u>i2 scales as the 
square root of the intersubband transition wave-length. 
In the ultra-strong coupling regime, the usually neglected 
anti-resonant terms of the light-matter coupling start 
playing an important role. In particular, the ground 
state of system is no longer the ordinary vacuum of pho- 
tons and electronic excitations, but rather a two-mode 
squeezed vacuum, whose properties strongly depend on 
the ratio $7^/^12. As this quantity can be dramatically 
tuned by applying an electrostatic bias, we have pointed 
out the possibility of observing interesting quantum elec- 
trodynamic effects reminiscent of the dynamical Casimir 
effect, i.e. the generation of correlated photon pairs out 
of the initial polariton vacuum state. A quantitative es- 
timate of the number of emitted photons has been given 
for the simplest case of an instantaneous switch-off of the 
light-matter coupling, and the results look promising in 
view of experimental observations. Work is actually in 
progress in the direction of extending the analysis to the 
case of a periodic modulation of flR.k, case in which one 
should be able to enhance the emitted intensity via para- 
metric resonance effects. From the theoretical point of 
view, this study requires a complete treatment of losses 
in order to describe the dynamical equilibrium between 
the parametric process generating the quantum radiation 
and the dissipation. 
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